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Abstract 

This article deals with a quantum-mechanical system which generalizes the 
ordinary isotropic harmonic oscillator system. We give the coefficients con- 
necting the polar and Cartesian bases for D = 2 and the coefficients connect- 
ing the Cartesian and cylindrical bases as well as the cylindrical and spherical 
bases for D = 3. These interbasis expansion coefficients are found to be an- 
alytic continuations to real values of their arguments of the Clebsch-Gordan 
coefficients for the group SU(2). For D = 2, the superintegrable character for 
the generalized oscillator system is investigated from the points of view of a 
quadratic invariance algebra. 



Paper to appear in the Proceedings of the VIII International Conference on 
Symmetry Methods in Physics (JINR, Dubna, Russia, 28 July - 2 August 
1997). The Proceedings of the Conference will be published in the Russian Journal 
of Nuclear Physics ( Yadernaya Fizika) . 



1 Introduction 



During the last 30 years, superintegrable dynamical systems have been the object 
of considerable interest (see [1-10] and references therein). In particular, numerous 
works have been devoted to the search for dynamical invariance algebras (especially 
quadratic algebras) of nonrelativistic systems with potentials presenting singular- 
ities. Such systems are important in various fields (e.g., Aharonov-Bohm effect, 
Dirac or Schwinger monopoles, confining problems, supersymmetry, etc.). 
It is the aim of this paper to investigate the system with the potential 

v=j:k, y^ = l^'-i+lp^^ p=ki-\ (1) 

a=l ^ ^ "''a ^ 

where i7 > and A;^ > (a = 1, 2, ■ ■ ■ , D). This system was already discussed for 
D = 2 by the late Professor Smorodinsky and his collaborators |]l| from a classical 
and quantum-mechanical point of view. We shall be concerned here mainly with 
D = 2 and 3 for which the spectrum of the Schrodinger equation 

Hm = E^>, H = -^A + V (2) 

shall be given. Emphasis shall be put on interbasis expansions in terms of analytic 
continuation of Clebsch-Gordan coefficients (CGc's) for the group SU(2). As another 
important result, we shall introduce a quadratic invariance algebra in the D = 2 
case. 



2 i^-dimensional case 

We briefly consider here the D-dimensional case in Cartesian coordinates. We start 
with D = 1 and look for a solution of the one-dimensional equation (^) for the 
potential Vi, see (||), with xi = x and ki = k. The resolution of this equation, with 
the conditions ^'(x) — > as x — > and oo, leads to the normalized wavefunction 



^n{x- ±k) = ^ ^^^'^ (v^)^" exp (-^x^j (fix^), neN (3) 

where L'^ is an associated Laguerre polynomial 0. The normalization is such that 

2 / ^Hn'ix; ±A;)*^„(x; ±k)dx = S^'n (4) 
Jo 

The discrete energy spectrum is given by 

E = Q{2n±k + l) 
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Only the sign + may be taken in front of k when k > ^. For < < |, both the 
signs + and — are admissible. For k = (|)~, due to the connecting formulas |TT[ 



between the (even and odd) Hermite polynomials 7ip{x) and the Laguerre polyno- 

lials Ln ^ {x^) and by f 
-, we immediately have 



mials Ln ^ {x^) and by putting p = 2n+ 1 for the sign + and p = 2n for the sign 



We now deal with the D-dimensional case. In this case, the Cartesian wavefunc- 
tion, that vanishes when — and oo (a = 1, 2, ■ ■ ■ , D), is 

D 



^„(x;k) = n <il^^{xa;±ka) 



a=l 

where n = ni ■ ■ ■ njj with Ua G N, x = Xi, ■ ■ ■ , x^) and k = ±ki, ■ ■ ■ , ^k^. The 
energy is 

r ^ 

E = n 2n + D + Y,i±ka) 

a=l 

where n = rii + n2 + ■ ■ ■ + nn is the principal quantum number. 

3 Two-dimensional case 
3.1 Cartesian basis 

In Cartesian coordinates ( = y), the wavefunction is 

±fci,±A;2) = ±A;i)^„2(y; ±A;2) (5) 

where \I'„^ (with a = 1,2) are given by (0). Note that we have the new constant of 
motion 



(in addition to the energy), where Da/s = —dap + ^ afi is the Demkov tensor |T2 . 



3.2 Polar basis 

In polar coordinates (p, v?), the potential (0) reads 

2 \ COS"' ip sm Lp 
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for which Eq. (0) may be separated by seeking a solution in the form R{p)^{ip). 
This leads to the system of coupled differential equations 



1.2 _ 1 

22 4 

cos^ if sin^ Lp 



^$ = 0, 



1 

dp (pdp) + 2E- ^fp^ - 



P 



R = (7) 



where ^ is a polar separation constant. 

The solution = *^',n('/'; ±^2) of the angular equation in (0) with the 

conditions 



is easily found to be 



(2m ±ki±k2 + l)m\T{m ±ki±k2 + l] 



X 



2T{m ± ki + l)T{m ± k2 + I) 
cos ^) 5^^=! (sin ^)i±fc2p^±fe2,±fci) (cog 2^) 



(9) 



where m G N and P^'^^ denotes a Jacobi polynomial. The normalization is such 
that 



Jo 

Then, the separation constant A is quantized as 

A = 2m±ki±k2 + l 
The radial solution R{p) = Rnpm{p] ±^1, ±^2) in (0) is 



(10) 



Rip) 



2nnp\ 



:i2) 



where rip G N is the radial quantum number. The function R satisfies the orthogo- 
nality relation 

/ ^n' m(P; ±^1. ±k2)Rn,m{P; ±^1> ±k2)pdp = 5^' 
Jo p P 

The energy corresponding to the n + 1 wavefunctions 

{p,<^]±ki,±k2) = R{p)H<^) 

(with n = Up + m fixed) is 

£; = ^](2n±A;l±A;2 + 2), nGN (13) 

where n is the principal quantum number. Note that only the sign + in front of 
ki and k2 has to be taken when ki > ^ and k2 > \. In the case < A;^ < | 



(with a = 1,2), Eq. (Qj shows that for each n we have four levels corresponding 
to (±^1, ±A;2). The degeneracy of the level with the principal quantum number n 
is n + 1. This degeneracy is identical to the one of the isotropic oscillator in two 
dimensions, for which the degeneracy group is SU(2). 

For ki = k2 = ily, we have A(i, i) = 2m + 2, A(-i, -i) = 2m and A(i, -i) = 
A(-i, i) = 2m + 1. Then, by using the connecting formulas |Tl| between Jacobi 
and Chebychev polynomials, we obtain the four following wavefunctions 



^2n,2m(P, V^) 
^I^2n+2,2m+2(p, V^) 



—^R2n,2m{p) cos2mip, U = 2n 



(14) 



'2n+2,2m+2(p) sin(2m + 2)(/3, n = 2ri + 2 (15) 



R 



'2n+l,2m+l(p) COs(2m + 1)(/?, Ti = 2n + 1 (16) 



2n+l,2m+l 



R2n+i,2m+i{p) sin(2m + l)v?, n = 2n+l (17) 



corresponding to the energy E = VL{n + 1). In Eqs. (|T^)- (|r^) , we have 

RpAp) = ^ (v^)*exp (-|/) Li-,(^]/) 

to be compared with the corresponding result for the ordinary circular oscillator. 
To close this section, let us mention that 



M 



1 
4 



^2 _ 1 

22 4 _^ 22 i 

cos^ 9? sin^ (p 



'/r2 _ 1 ^2 _ 1 
22_ 4 _|_ 22. i 



is a polar constant of motion, the eigenvalues of which are A (see ([TTl) ). 



(18) 



3.3 Connecting Cartesian and polar bases 

According to first principles in quantum mechanics, we have 

n 

*n,n. = E W^n"^n.(±fcl> ±h)^n,m (19) 
m=0 

where Up + m = ni + n2 = n. In Eq. (|19D, it is understood that the wavefunctions 
both in the left- and right-hand sides are written in polar coordinates {p, tp). Fur- 
thermore, by using the asymptotic formula for the associated Laguerre polynomials, 
Eq. (|l^) yields an equation that depends only on the variable p. Thus, by using the 
orthonormality property of the function $ with respect to the quantum number m, 
we obtain 

WZn,{±K±k2) = {-irBZn,i±h,±k2)EZn,{±h,±k2) (20) 
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where 



and 



f22) 



[n 



m)!m!r(m ± A;i ± A;2 + l)r(n + m ± A;i ± A;2 + 2) 



\ ni!n2!r(m ± fci + l)r(m ± A;2 + l)r(ni ±ki + l)T{n2 ±^2 + 1) 



By making the change of variable x = cos 2ip and by using the Rodrigues formula |Tl 
for the Jacobi polynomial, Eqs. (pOD-(p2D lead to the integral representation 



(2m ± A;i ± A;2 + l)(n - m)\T{m ±ki±k2 + l)T{n + m±ki±k2 + 2) 
\\ niWmlTim ± ki + l)T{m ± k2 + 2)r(ni ± ki + l)T{n2 ± ^2 + 1) 



X 



2ni+n2+m±k-i±k2+l J _^ 



1 



X 



I "2/ 



1 



X 



dx'^ 



-[{l-x 



1 + X 



,m±A;i 



(23) 



for the interbasis expansion coefficients WJ^^j^^{±ki,±k2). 

Equation (^) can be compared with the integral representation [|T3| for the 
CGc's {aba(3\c'j) of the group SU(2). This yields 



M^-„^(±fci,±fc2) 



4)"^-"''(a6a/?|c7) 



(24) 



with 2a = ni + 722 ± ki, 26 = tt-i + 77-2 ± A;2; 2c = 2m ± /ci ± ^25 2q! = rii — 77-2 ± A;i 
and 2(3 = n2 — k2- Since the quantum numbers in (|2^ ) are not necessarily 
integers or half of odd integers, the coefficients for the expansion of the Cartesian 
basis in terms of the polar basis may be considered as analytical continuation of the 
SU(2) CGc's. 

The inverse of Eq. (plQ]), viz.. 



y: w:u±k,,±k2)^ 



nin2 



(25) 



ni=0 



follows from the orthonormality property of the SU(2) CGc's. Thus, the relation 

W:U±ku±k2) = WZn2{±k,,±k2) 



gives the expansion coefficients in (g^). The SU(2) CGc's can be expressed |T3 

WZn2{±h,±k2) = {-ir 



m 



terms of the hypergeometric function 3^2(1), so that Eq. ( P^ can be rewritten as 

?T,!r(?T,2 + m ± ^2 + 1) 



npWim ±ki + l)r(m ±^2 + 1' 



(2m ± A;i ± A;2 + 1) 



T{ni ±ki + l)r(m ±ki±k2 + l) 
ni\n2\m\r {n2 ± k2 + l)r(n + m ± fci ± ^2 + 2) 

n — m^ki^k2 — 1, —7^25 —m 



\ 



-ni - n2, - n2 -mT k2 
By using symmetry properties for 3F2(1), we arrive at the expression 

(-l)'^n! 



w::^(±k,,±k2) 



r(i ± k2) 



\ 



(2m ±ki±k2 + l) 



r(m ±ki±k2 + l)r(m ±k2 + l) 
ni\n2\m\np\r{m ± fci + 1) 



X 



T{ni ±ki + l)r(n2 ±k2 + l) 



\ V{n + m±ki±k2 + 2) ^ '^\l±k2, -n^-n2 



-m, m ± ki ± k2 + I, — n2 



(26) 



Alternatively, by using the formula [T^ connecting the Hahn polynomial h^^'^^ and 
the function 3^2(1), we obtain 



wz,J±h,±k^) = (-1)' 



(2m ± A;i ± A;2 + 1) 



m!np!r(m ± A;i ± A;2 + 1) 
ni\n2\T[m ± fci + l)r(m ± A;2 + 1) 



r(n,±fci + l)r(n2±fc2 + l) ni+n. + l) 

^ r(n + m±A;i±A;2 + 2) ^"^ [n,,n,+n2 + l) 

in terms of Hahn polynomials. 
3.4 Invariance algebra 

Let us consider the following realization of the SU(1,1) generators 



(a) 



2 1 



_ 

2^2 I a 4 



2 7(1) 



We thus have two copies (for a = 1, 2) of the Lie algebra SU(1,1) given by 

r jia) jia)-\ _ ■ j(a) r j(a) j(a)-, _ ■ j{a) r j{a) j{a)i _ ■ j{a) 
["^0 1 "^l J ~ '"^2 ) ["^1 ) "^2 J ~ '"^0 ) ["^2 ) "^0 J ~ ""^l 

with the Casimir operator 

Qa = [Jt^r-[4''^r-[jt^r = \{ki-i) 



(27) 
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Introducing the raising and lowering operators J^*^ = j'f'^ ± 132^., we get 

J^)] = ± [Ji"\ Jf^] = 2 J^'^) and = J^^ji"^ 

As an irreducible representation of SU(1,1), the positive discrete series consists of 
an infinite number of states. Each of these states will be denoted as \ jaTT^a), where 
"la = ja + na {jia = 0, 1, 2, ■ ■ ■). The eigenvalue of the Casimir operator is 

Qa = jaija " 1) 

SO that from ( P?] ) we have ja = |(1 ± ka)- The matrix elements of the generators of 
the group SU(1,1) may be obtained through 



J^^jama) = rUaljama), J^^\jama) = (ma±ia) (m-a=Fia±l) |ia?^a±l) (28) 

with J^'^ljaja) = 0. Let us now define 

Co = 4'' + 4'\ C^ = 4' + jS' (29) 



Equation ([29| ) corresponds to the direct sum of the two SU(1,1) algebras for a = 1, 2. 
The coupled basis \jm) satisfies 

Co\jm) = m\jm) = {j + n)\jm), Q\jm) = j{j - l)\jm) 

Given the values ji and j2, the parameter j can take the discrete values 

j = ji + J2 + g, g e N 

The Clebsch-Gordan decomposition yields 

\jm) = J2 {jij2mim2\jm) \jimi) (g) \ j2m2), m = rrii + 1112 

mim2 

with 2ja = I ± ka, Inia = 2na + 1 ± A;^ and 2j = 2q + 2 ± ki ± k2- By using the 
connection between the SU(1,1) CGc and the 3^2(1) function |15[, one can obtain 
the same hypergeometric function as in (|26|). 

Note that the Hamiltonian H of our two-dimensional oscillator system is 

H = 2nCo 

From ( PBD and (PUj), we recover the spectrum of the system as given by (|T^ with 
n = rii + n2. 

Let us consider the two following operators 

N = - 4'\ M = Q, + Q2 + 24^) - j« ji:^) - Ji^) Jf + \ 
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They commute with H. Indeed, they are nothing but the integrals of motion 
and (|T^). Moreover, let us define a third operator T via T = [N, M]. We have 



T = 2 



or 



T 



1 



ydy + 



1\ k. 



kf - - 



2 _ 1 

2 4 



a;9:r + 



The operators A^, M, T and span a closed quadratic algebra since 
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[M,T] = -2(MA^+ArM) + '"^ ^^'^ H-N, [T,N] = -2N'^ + ^H'^ -AM-kl-kl-l 

hold in addition to [iV, M] = T, [A^, i/] = and [M, //] = 0. 

In the limiting case /ci = A;2 = |, we obtain a quadratic algebra too. In this case 

Instead of A^, Li and T, we can consider A^, and [A^, L^] . In this regard, by 
putting 

Pi = N, P, = P, = = \^D^y 

we end up with the Lie algebra corresponding to the commutation relations 

[Pk, Pt] = ieumPm, A;,£,m G {1,2,3} 
Finally, going back to the generic case for ki and k2i we define 

Lo = ~{2NTki±k2) 



and 



'{ni ± ki){n2 ± /C2 + 1) ^in2 ± A;2)(ni ± A;i + 1 

They act on the eigenfunctions (|^) of the Hamiltonian H as 

1 ' 



1 i\ / 1 r, 

^1 + 2 T 2 j ^ 2 ^ 2 )*"iTi"2±i 



The operators Lq, L+ and L_ generate the Lie algebra SU(2) with 

[Lo,L±] = ±L±, [L+,L_] =2Lo 
and are closely connected to our integrals of motion. 
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4 Three-dimensional case 



4.1 Spherical basis 

In spherical coordinates (r, 9, ip), the potential (|l]) can be rewritten as 

1 



2 2r2 



/ iL2 _ 1 1,2 _ 1 \ t.2 _ 1 
21 4 _^ ^^2 4 \ _^ 23 4 

sin^ 5 I cos^ (p sin^ 99 / cos^ 9 



Looking for a solution of Eq. (|^) in the form R(r)Q{9)^{ip), we are left with the 

system 



sin 9 



de{sm9de) + J{J + r 



h.2 1 
'^1 4 


1.2 1 
'^2 4 


cos^ 


sin^ (p 




h.2 1 
'^S 4 


sm^9 


cos^ 9 


- Q^r' - 


J(J+1) 





e 



R 







(30) 
(31) 
(32) 



The solution $(95) = $m(v^; ±^2) of Eq. (pOD, satisfying the boundary con- 
ditions (P) and the normalization condition (p!OD, is given by (P). The separation 
constant A in ( |5(]| ) and ( pi] ) is quantized according to (0). 

The solution Q{9) = Bg^niO; ±ki ±k2± h) of (P) is (see [|) 



[2(m + g + 1) ± A;i ± A;2 ± k3]q\T{q + 2m ± ki ± k2 ± + 2) 



T{q ±k3 + l)T{q + 2m + 2±h±k2) 

X (cos 9)'^^''' (sin efpi^'^^s) (cos 2^) 
which satisfies the boundary condition 



0(o) = e(-) = o 

and the normalization condition 

2 e,,™(^; ±A;i, ±A;2, ±A;3)*e,„(^; ±A;i, ±A;2, ±^3) sin^d^ = 6, 

'J 



The spherical separation constant J in (^) and (|32D is 



1 3 

J = 2q + A±k^ + - = 2q + 2m±ki±k2±hi + - 



The solution R{r) = Rn^qmij] ±A;i, ±^2, ±A;3) of Eq. ( ^2]) is 
i?(r) 



2^271^! 



^ r(n. + 2g + 2m±fc,±fc2±fc3 + 3) V^)' ("i'^') 



with 



/•OO 



where ri^ € N is the radial quantum number. 
The energy of the system is 

3 

E = n{2nr + J+-) = n{2n ± ki ± k2 ± k^ + 3), n G N 

where n = Ur + q + m is the principal quantum number. It corresponds to the 
wavefunctions 

^nrgmir, 0, ±k,, ±k2, ±k3) = R{r)Q{e)^ip) 

with n fixed. 

4.2 Cylindrical basis 

In cylindrical coordinates {p,ip,z), we have 



2 ^ +2p2l^cosV sinVj + 2r^ + 



The corresponding Schrodinger equation may be solved by looking for a solution in 
the form R{p)^{Lp)Z{z). By combining the results of Sections 2 and 3, we get 

Z{z) ^'^nM^h), <l>i^) = <l>U^;±ki,±k2), R{p) = RnpUp;±h,±k2) 

as given by (H), (H) and (0), respectively. The energy 

E = n{2n ±ki±k2±k3 + 3) 

corresponds to the wavefunctions 

(p, If, z; ±ki,±k2, ifcs) = R{p)^{if)Z{z) 

for which the principal quantum number n = Up + m + is fixed. 

4.3 Connecting Cartesian, cylindrical and spherical bases 

In the three-dimensional case, we have 

rii+ri2 np+ns 

^nmana = ^J^na (=^^1' =^^2)^ J2 ^npnsi^ki,±k2,±h)'^nrqm 

m=0 q=0 

where rii + n2 = m + Hp and rir + q = Up + n^. For the expansion of the Cartesian 
basis over the spherical basis, we have 

^n.n.na = E CZl.n,i±h, ±k2, ±ks)^n.,m (33) 
mq 
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where ni + 712 + % = Ur + q + m. The coefficient ^^^^^(iA;!, ±^2) is identical to 
the one found in the two-dimensional case. It is given by Similarly, it is easy 
to obtain 

Vl^.^i±h,±k2,±h) = (-l)"''-«(a'6V/?'|cY) (34) 

where 2a' = na + rip ± k^, 2b' = + Up + 2m + 1 ± /ci ± A;2, 2c' = 2q + 2m + 
1 ± /ci ± A;2 ± /cs, 2a' = — Up ± and 2(3' = 2m + Up — + 1 ± ki ± k2- The 
expansion coefficients in (|33D are given by the formula 

CZl,nMk,,±k,,±k,) = w:U±k,,±k,)Vl^.^^i±k,,±k,,±k,) (35) 

The value of the right-hand side of (pSj) follows from (^) and (34). 
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